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Density waves in the shearing sheet 
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Abstract. The shearing sheet model of a galactic disk is studied anew. The theoretical description of its dynamics is 
based on three building blocks: Stellar orbits, which are described here in epicyclic approximation, the collisionless 
Boltzmann equation determining the distribution function of stars in phase space, and the Poisson equation in 
order to take account of the self-gravity of the disk. Using these tools I develop a new formalism to describe 
perturbations of the shearing sheet. Applying this to the unbounded shearing sheet model I demonstrate again 
how the disturbances of the disk evolve always into 'swing amplified' density waves, i.e. spiral-arm like, shearing 
density enhancements, which grow and decay while the wave crests swing by from leading to trailing orientation. 
Several examples are given how such 'swing amplification' events are incited in the shearing sheet. 
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1. Introduction 

The shearing sheet (Goldreich & Lynden-Bell 1965, Julian 
& Toomre 1966) model has been developed as a tool to 
study the dynamics of galactic disks and is particularly 
well suited to describe theoretically the dynamical mech- 
anisms responsible for the formation of spiral arms. For 
the sake of simplicity the model describes only the dynam- 
ics of a patch of a galactic disk. It is assumed to be in- 
finitesimally thin and its radial size is assumed to be much 
smaller than the disk. Polar coordinates can be therefore 
rectified to pseudo-cartesian coordinates and the velocity 
field of the differential rotation of the disk can be approx- 
imated by a linear shear flow. These simplifications allow 
an analytical treatment of the problem, which helps to 
clarify the underlying physical processes operating in the 
disk. Goldreich & Lynden-Bell (1965) used hydrodynam- 
ical equations as dynamical equations, which were com- 
bined with the Poisson equation in order to take account 
of the self-gravity of the disk. Their principal result was a 
gravitational instability giving rise to spiral-arm like den- 
sity enhancements, which take part in the general shear 
of the disk and grow while swinging by from leading to 
trailing orientation. This effect has been pursued and dis- 
cussed further in many studies, notably by Goldreich & 
Tremaine (1978). Independently Julian & Toomre (1966) 
have investigated a shearing sheet made of stars by us- 
ing the collisionless Boltzmann equation combined with 
the Poisson equation and found also 'swing amplification' 



of spiral density waves with the additional effect that the 
density perturbations die out eventually due to smear- 
ing out of the ever tightening spiral arms by the finite 
epicyclic motions of the stars. Toomre (1981) has argued 
vigorously that this swing-amplification mechanism is the 
principal process of spiral arm formation in galactic disks 
and that the spiral arms of many spiral galaxies are thus 
of transient nature. Indeed, numerical simulations of the 
dynamical evolution of self gravitating, differentially ro- 
tating disks by Sellwood & Carlberg (1984) have shown 
impressively recurring swing-amplification events, which 
lead to an ever-changing appearance of the disk. 

The alternative concept of quasi-stationary density 
waves has been pursued in the very influential work of Lin, 
Shu, Bertin and collaborators culminating in the papers 
by Bertin et al. (1989 a,b) and Lowe et al. (1994). In their 
work they always consider full disk models. One key ingre- 
dient of their recent models is the introduction of an inner 
reflecting Q-barrier, which shields the inner Lindblad res- 
onance. Together with the semi-reflecting corotation zone 
this forms a kind of resonance cavity, in which the rigidly 
rotating density waves grow exponentially. 

A general theory of modes of star disks has been devel- 
oped by Kalnajs in a series of papers (Kalnajs 1971, 1977, 
1978), which has been used by him and several other au- 
thors, notably Hunter (1992) and Pichon & Cannon (1997) 
to describe spiral and bar-like perturbations of the disks. 
Density waves as modes of disks have been also searched 
for and localized in numerical simulations of the dynami- 
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cal evolution of star disks by Donner & Thomasson (1994) 
and Scllwood and collaborators (see Sellwood & Evans 
2000 and references therein, also Toomre 1981). Toomre 
and Kalnajs (1980) have introduced inner boundaries into 
their numerical, softened gravity, simulations of the evolu- 
tion of the shearing sheet and describe two kind of modes 
growing in the sheet. The present paper is the first of 
a series meant as a complement to these studies. I con- 
sider again a shearing sheet made of stars and discuss 
first the unbounded sheet and then the dynamical con- 
sequences, when inner boundary conditions are applied. 
The aim is to give a consistent theoretical description of 
pure swing-amplification as well as exponentially growing 
modes in the framework of the same model. In this paper 
I consider the unbounded shearing sheet. The purpose is 
to develop the formalism, and I obtain mainly previously 
known results again. Contrary to previous studies I avoid 
Lagrangian - shearing - coordinates and stay in the Euler 
picture. This leads to a compact, transparent form of the 
central Volterra integral equation (cf. section 6), and will 
enable me also later, when I consider a bounded shearing 
sheet, to derive a dispersion relation for rigidly rotating, 
exponentially growing modes of the disk. 

2. Stellar orbits 

The first step required in order to construct a stellardy- 
namical model of a galactic disk is to describe the stellar 
orbits in the disk. The equations of motion of the stars 
are derived from the Lagrangian 



(1) 



where r, 8 denote polar coordinates. Since the disk is as- 
sumed to be infinitcsimally thin, only plane orbits are con- 
sidered, r, 8 are the corresponding radial and angular ve- 
locities of the stars. <&(r, 8, z) is the gravitational potential 
of the disk. The potential is assumed to be axisymmetric, 
$(r, z), in the present section. Non-axisymmetric pertur- 
bations of this basic state of the disk will be taken into 
account in the next section. Considering now a patch of 
the disk as in the shearing sheet model I introduce pseudo- 
cartesian coordinates 

x = r — r and y — r (8 — Clot) . (2) 
ro denotes the galactocentric distance of the center of the 



patch and Qq 



r V dr ) r 



is the mean angular veloc- 



ity of the patch around the galactic center. In order to 
keep the analysis tractable I consider the stellar orbits 
in epicyclic approximation (Lindblad 1959). This assumes 
that 



\x\ <C r and \y\ -C fio^o • 



(3) 



Taylor expansion of the Lagrangian with respect to x leads 
to 



1 
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y 



4fl xy + AACl x 2 ] , (4) 



where A denotes Oort's constant, A = — ^ro (^) r ■ 
The equations of motion are then given by the standard 
epicyclic equations, 



x = 2fl y + AACIqx and y = —2CIqx , 



(5) 



which can be solved straightforward in terms of linear 
and trigonometric functions (cf. equations 10 and 11 be- 
low). The linear terms give the mean motion of the stars, 
whereas the trigonometric terms describe the epicyclic os- 
cillations of the stars around the mean motion with a pe- 
riod of k = ^f—AVl^B, where B = A — fl denotes Oort's 
second constant. The plane stellar orbits conserve two in- 
tegrals of motion, one related to angular momentum, 



dy 



Cl r + 2fl x + y, 



and the radial action integral, 



4B 2 



J1 = 2- K {U 



(6) 



(7) 



The radial action integral is defined in terms of the radial 
velocity of a star u — x and the circumferential velocity 
referred to the mean circular velocity of all stars at the 
position of the star, v = y + 2Ax = J 2 — Cl a r a + 2Bx. Next, 
the Hamiltonian of the stellar orbits is derived from the 
Lagrangian (4) and expressend in terms of the integrals of 
motion giving 



H = KJl + 2l? (J2 ~ °° ro)2 ~ l^ 



(8) 



Angle variables w\, W2 are introduced as canonical conju- 
gate variables to the integrals of motion, 



Wi 



dH 
8J[ 



K , W 2 



dH _ A 

dI>~B 



(J 2 - fl r ) . 



(9) 



Assembling the results from above, the orbit equations can 
be written as 



J 2 - Slo^o 2Ji . 
x = —^- + ] / — smWl , 



y = w 2 



^f2lUl 



COSW! , 



2B 

u = COS Wi , 

2B /77—r ■ 
v = — vzkJi sinwi . 

K 



(10) 



These equations show that the mean motions of the stars 
are given by 



J 2 - ft ro 

x = — — = const. , 

— 2B 

A 

y = —(J 2 -Clor )t = -2Axt, 



(11) 



describing the linear shear flow. The trigonometric terms 
describe the epicyclic oscillations around the mean mo- 
tion, with their amplitude determined by the radial action 
integral J\. 
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3. Boltzmann equation 

The ensemble of stars is described statistically by the dis- 
tribution function of stars in phase space /. Its time evolu- 
tion is determined by the collisionlcss Boltzmann equation 



df 
dt 



[/,#] = 0, 



(12) 



where the Poisson bracket indicates the divergence of the 
Hamilton flow in phase space. If integrals of motion are 
used the Poisson bracket has the form 



[f,H] 



df dH df dH 



dw\ dJi dw2 9 J 2 



(13) 



Following Julian & Toomre (1966), the Schwarzschild dis- 
tribution, 



fo 



S 



-2B 27177,2 



exp - 



kJi 



(14) 



is chosen as the distribution function of the homogenous 
background of the sheet. E denotes the surface density of 
the disk, which is kept constant over the sheet. a u is the 
radial velocity dispersion of the stars. The circumferential 
velocity dispersion is given according to the epicyclic ra- 
tio <t v /<t u = AB 2 /k 2 . Since the stationary Schwarzschild 
distribution depends only on the radial action integral, it 
solves the Boltzmann equation by definition (cf. equation 
13). 

4. Linear perturbation analysis 

Spiral arms are usually thought to be only minor pertur- 
bations of galactic disks. Thus I choose a perturbation 
'Ansatz' of the form 



f = f + Sf, H = H Q + 6<1>. 



(15) 



where H is the Hamiltonian (8) of the unpertubed or- 
bits. The physical idea is to subject the shearing sheet to 
a small spiral-like potential perturbation 5$ and to de- 
termine the disk response 5f by solving the Boltzmann 
equation. The Boltzmann equation is linearized accord- 
ingly, 



+ [fo,S*] + [5f,H Q ] = 0, 



(16) 



dt 

where the quadratic Poisson bracket [5f, (5$] has been ne- 
glected. The role of the quadratic term will be investigated 
in a further paper (Fuchs 2001, in preparation). A brief 
first report of that work was given in Fuchs (1991). The 
potential perturbation is Fourier transformed with respect 
to the spatial coordinates and I consider a single Fourier 
component 



(5$ k = $ k exp i (k K x + k y y) 



(17) 



in equation (16). The Poisson brackets have then the ex- 
plicit form 



and 



fl r ] 



dw 2 ' 



(19) 



which leads to the linearized Boltzmann equation 



dSf 86 f A. T 



o ^ 95 f 

OW 2 



+ 



ik x cos toi + — ikysinwi 
2B 



So 
2-Kat 



■ exp 



) $ k exp (i/c x x 

K ) 



2B 
ik y y) = . 



(20) 



When x and y are substituted from equations (10), equa- 
tion (20) is a linear inhomogenous partial differential equa- 
tion with t, wi, and W2 as independent variables. It be- 
comes immediately clear that the coefficients of the equa- 
tion do not depend on the time coordinate t and the angle 
variable W2 ■ This allows a Fourier transformation of equa- 
tion (20) with respect to t and w 2 , which may be written 
in abbreviated form as 

iv8fu,i + k ^/^- + 4 {'h - ^or ) U5fu,i 
dw\ B 



/oo 
dw' 2 e- ilw ^Inhom u (Ji, J 2 , wi) e 
-OO 

= 2ne aw2 Inhom^ (J 1; J 2 , Wl ) 6 (k y - I) 



(21) 



where S(k y — I) denotes a delta function. Equation (20) 
will be solved explicitely below. However, it becomes clear 
from equation (21) that the solution must be propor- 
tional to &ke tlW2 8 (k y — I), because of the linearity of the 
equation and the consideration that only a non-zero per- 
turbation leads to a disk response. Transforming back to 
u> 2 space, this implies a functional dependence of Sf on 
u> 2 as 

Sfoce ik » W2 . (22) 

Furthermore it is advantageous to split off the background 
distribution function from the disk response, 8f u = fi lOJ fo, 
and I obtain the ordinary differential equation 

4fi,u 



V2^7T 
■exp{ik x 



dw\ 
ik x cos wi 

J 2 - ^07Y) 

-2B 



A , 

iky— (J 2 



^Vo) fl,u 



2B 



iky sin w 



+ 



2Ji 



sma)! 



y 2B C0SWl J 



0. 



(23) 



Equation (23) can be solved by 'variation of the constant'. 
Very similar equations have been treated with the same 
method by Lin, Yuan, & Shu (1969) and Julian & Toomre 
(1966) and I refer in particular to the lecture notes by 
Lin (1970) for a full explanation of the technique. The 
homogenous part of equation (23) has solutions 



ik x cos wi 



2B 



iky sinwi 



./o<5$ k (18) /i^ocexp- 



w + ^k y (J 2 - Qor ) ) wi 



(24) 
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and 'variation of the constant' means to look for a partic- 
ular solution of equation (23) in the form 

fi,u> = C(wi) ex P-^ + Jj k y i'h - ^o^o)^ Wi . (25) 

Inserting this into equation (23) and solving for C(w\) 
gives 

/ \ i£, ( ., Ji-^Tq\ ._ 
C{wi) = -exp z/c x — — )^ kiU1 expir]w 



J 

Jo 



Wl — w 



-2B 



dw[ cos w[ exp i [£ sin w\ + r)w'^\ , (26) 



where the abbreviations £ = y k 2 + jj^k 2 and 77 = 

i (« + fc y | ( J2 — fio r o)) have been introduced to keep the 

formulae compact, w is defined asS= arctan^^^^. 

The general solution of equation (23) is then given by 
the particular integral (25) plus the solution of the ho- 
mogenous equation multiplied by an arbitrary integration 
constant, 



<fk. 



1,LU 



exp ik x 



-2B 



{D — i£ exp i-qw (27) 



f 

Jo 



dw[ cos w[ exp i [£ sin + rjw^}} exp —irjw\ . 



The integration constant D is determined by the require- 
ment that the distribution function must be defined in a 
unique way in velocity space, i.e. must be periodic with 
respect to the angle variable W\ with a period of 2tt, 



fi,u (wi) = fi,u> (wi + 2ir) 



(28) 



Some algebra as indicated in Lin (1970) leads to the final 
result 



fi 



-^exp (ik x — — jyp?- )<;»!, , < sp(/'i sin (ir t - ir)) 



v f +7T 

■ 1 1 ; — - / dw\ expz \riw\ — £ sin w\ 

X 2sin(7rr,) J _^ 1 

• cos (w\ — w) — £ (1 + cos w^) sin (w\ — w)] } . 



(29) 



The disk response to the single Fourier component of the 
potential perturbation is thus given by 



Sf = /o/i,w exp (iwt + ik y w 2 ) 



(30) 



Note that according to equations (10) and the definitions 
of £ and w 

kx— — ^° r ° + £, sin (wi —w) + k y w 2 = k x x + k y y . (31) 
— 2B 

This shows that the disk response follows to a first approx- 
imation the spatial pattern of the potential perturbation. 
However, the second term in the curly brackets in equa- 
tion (29) depends through r\ on (J 2 — ^o?~o) an d thus on 
the spatial coordinate x, so that contrary to <&k,w, fi,u is 
not a Fourier component of a Fourier decomposition of 
the disk response. 



5. Self consistent perturbations 

The shearing sheet is assumed to be self-gravitating. Thus 
the potential perturbation and the corresponding disk re- 
sponse have to satisfy the Poisson equation, 



{ 



d 
dx 2 



d 
dy 2 



d 
dz 2 



}<5$ = AirGSZS (z) , 



(32) 



where <5E denotes the surface density of the disk response 
and 5(z) is a delta function with respect to the verti- 
cal spatial coordinate z. G is the constant of gravity. 
Background potentials do not enter into equation (32) be- 
cause of the linearity of the Poisson equation. Off the mid- 
plane the right hand side of equation (32) is equal to zero. 
This leads immediately to a three-dimensional extension 
of a potential perturbation (17) in the form 



<5$k = $k exp i (k x x + k y y) - ^Jk^ + k 2 



(33) 



Integrating both sides of the Poisson equation with re- 
spect to z gives the surface density required to sustain the 
potential (33), 



+00 02 Q2 Q2 



' dx 2 dy 2 dz 2 J 



■ exp 



i (k x x + k y y) - Jk 2 + k 2 \z\ 



= -2y k 2 + k 2 $ k exp i (k K x + k y y) 
= 47rGSk exp i (k x x + k y y) , 



(34) 



where £k denotes the Fourier coefficients of the Fourier 
transform of ST,. 

It was pointed out in the previous section that the disk 
response to a single Fourier component of a general po- 
tential perturbation as considered in equations (17) and 
(33) is not a Fourier component of the disk response to the 
general potential perturbation. Kalnajs (1971) has shown 
a way to overcome this difficulty by taking Fourier trans- 
forms of both sides of the Poisson equation in the following 
way 



1 



(2nY 



/+00 r+00 r+00 
dz dx dye' l « x+k 'y y ) 

-OO J —CO J — OO 



d 2 



^ dx 2 ^ dy 2 dz 2 



/>+oc r+00 
dk x / , 
-00 J —00 



• exp 



i (k x x + k y y) 

+oc z' + oo 



t'2 

h v 



k 2 \z\ 



= il^l f dx f dye-^+Kv) 
(2ir) J -co J -co 



(27T) 

dk x / dk y / du / dvfoh,„e ik » w * . (35) 

-00 J -co J —co J —co 

The k x ,k y integrations indicate that all Fourier coeffi- 
cients have to be considered together, while the integra- 
tions of the distribution function of the disk response over 
the u und v velocity components on the right hand side of 
euation (35) give the surface density of the disk response. 
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The left hand side of equation (35) can be evaluated by w\ integration of the first part of the integrand leads to a 
carrying out first the quadrature with respect to z, giving quadrature 



+00 



dx 



+00 



dy 



+00 



dk* 



+00 



dk v 



k 2 + k 2 $ k exp i [(k x -k' x )x+ (fc y - k' y ) y] 



f 2 ' 
Jo 

-I 



dwi exp i [£ sin (wi — w) — £' sin (wi — w')] 



= -2^k' 2 + fc y 2 $ k , 



(36) 



d Wl expi[(Z 2 - 2^'cos(w -W) + £ ' 

• sin (wi — 1 



(40) 



Equation (35) represents in effect an integral equation 
which has to be solved to obtain self consistent pairs of 
potential and density perturbations of the sheet. 

6. Derivation of the Volterra integral equation 



where £' and w' are defined analogous to £ and w for 

primed wave numbers, w is given by 

w = arctan ((£ sin w — £' smw') / (£ cos w — £' cos zzl')). 

A Fourier expansion of the integrand of (40) with re- 
spect to w\ (Abramowitz & Stcgun 1968, section 9.1 ) 
gives immediately the result of the quadrature 



In order to derive the central integral equation of this 
study the right hand side of equation (35) has to be evalu- 
ated explicitly. The multiple integrals imply an integration 
over 

the integrations from spatial x, y and u, v velocity coordi- 
nates to Ji, J 2 and wi, w 2 variables. Using equations (10) 
it can be shown that the Jacobian of the transformation 
is equal to one, 



(41) 



where Jo denotes the Bessel function of the first kind, 
the entire P has7space"and\°UsVonvcnient"to change The second P art of the int egrand of (38) is slightly more 



complicated, 

,-2jv 



/ dw\ exp i [— £ sin (w[ + w\ — w) — £' sin (wi — w')] 
Jo 



d(x,y,u,v) 



d(Ji,J 2 ,wi,w 2 ) 



2tt J ( (£ 2 + 2££' cos (w[ -w + Vf) + £ 



/2\ 2 



1. 



(37) 



(42) 



I consider next the integration over the action variable Ji 
which leads for the first part of the integrand of (38) to a 



Inserting / , fi, u , x, and y into the right hand side of quadrature 
equation (35) leads to 

d Ji exp - 



4irG 



(2tt)' 
So 



2tt 



dw\ 



+ 00 



2na 2 
■ expi 

•{1- 



exp - 



K T 



dw 2 
+00 



dJ 2 



r+°° $ k 
i dky — 2" 

J — 00 



Jl 2n J 



K 1,2 
40 2 y 



-2B 



-\lk 2 + 



K yi 

ab 2 y 



, J2 - £Vo , . t . , _^ 
^ — 1 - k y w 2 + £ sin (wi — w) 



-2B 



n 



+00 



2 sin (7777) 
cos (wi — w) — £ (1 + cos w[) sin (w\ 



dw' 1 expi[r)w' 1 — £ sin w[ 

m 



■\l k 'x + ^2 k y cos & - ™') ] 5 V^) . (43) 

if £ and £' are inserted into the argument of the Bessel 
function. According to formula (6.631) of Gradshteyn & 
Ryzhik (1965) this can be evaluated as 



2na 2 



•exp — 1 



, J 2 - Qpro 
x -2B 



-exp 



1 



k' x \ sinwi 

K 



k 2 4- — k 2 

x 40 2 y 



+ 



1 



k' 2 + — k' 2 

x 40 2 y 



, , V2kJi 

+fc y w 2 - fc y 2 ^ coswi 



(38) 



fc2 



— k 2 

4B 2 y 



— fc' y cos («;-«/) 



(44) 



and similarly for the second part of the integrand of (38) 



All of the quadratures but one are, even if tedious, 
straightforward. The integration over the angle variable 
w 2 is particularly simple, 



2k<t 2 



-exp 



1 

2^ 



+ OO 



dw 2 e< k *- k 'y) W2 =5(k y - fc y ) , 



(39) 



ui 1 1,2 

x 402 K y 



fr 2 4- — — k 2 

x 402 y 



k' 2 4- — — /!r' 2 
402 y 



A;' 



^fc' 2 
402 K y 



• COS (w[ -10 + uf) 



(45) 



and the integration over k y contracts the integrand of (38) 
to fcy = fcy. This shows that the integral equation (35) is 
separating in the circumferential wave numbers fc y . 

The integrand of equation (38) is split according to 
the two terms in the curly bracket into two parts. The 



In the second part of the integrand of (38) a quadrature 
over the auxiliary angle variable w[ has to be carried out, 



+7T 



dw[ expi [r)w[ — A cos (w[ — w + w')] , 



(46) 
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wher e I have in tr oduced A as an abbreviation for A = 
\Jk%. + -^k y ^Jk' x + j^j-fc'y". A Fourier transform with 
respect to w[ gives (Abramowitz & Stegun, 1968, Section 
9.6) 



/+7T 



dw[ (cos^Wi) + isin^wi)) 



Io (-A) + 2 ^ (-A) cos(n (w[ - w + w')) 



2 sin(7T7?) 



+2 Y In (-A) (-1)" sin (nrj) cos (n (u; - w')) 



n=l 
1 1 



rj — n rj + n 

oo 

+2i J n (-A) (-1)" sin (7177) sin (n (W - w')) 

71=1 
1 1 



?7 — n 77 + n 



(47) 



where the /„ denote modified Bessel functions of the first 
kind. 

This leaves finally the integration over the variable J2 . 
The first part of the integrand of (38) is evaluated easily, 



/+00 
- 00 



d Ji exp 



i (k x k x ) 



J 2 - fVo 



2tt (-2B) exp 



i (k x k x ) 



-2B 
25 



<5(fc x -fc x ). (48) 



The integration over fc x contracts the first part of the in- 
tegrand of (38) to fc x = k' x . 

However the second part of the integrand needs more 
care. The first term of the series (47) leads to 



/+00 
-00 



d J2 exp 



i [k x k x ) 



J 2 - Hon) 



= -2tt (-2B) exp 



-2B 

i {k-x. k x ) 



2B 



(-l)/o(A) (49) 
Io (A) 6 (k x - k'J , 



and the integration over k x contracts this part of the sec- 
ond half of the integrand again to k x ~k' x . The next terms 
of the series (47) lead to integrals of the form 



/+00 
-OO 



■(-»?) 



d J2 exp 
1 



l {k x -K)' h -^ 



-2B 



± 



1 



rj — n rj + n 



(50) 



The auxiliary variable rj is substituted for the variable J2 

and we have drj = ^-x-dJi- Since the Oort constant B is 
negative, I consider first the case of negative wave numbers 
k' y . The integral (50) is of the type 



lim 

A^O 



1-00+2 



-00+- 



drje 



iAr; V_ 



77 ± n ' 



(51) 



where A denotes A = k (k' x — k x ) /(2Ak' y ). According to 
its definition the variable i] is related to the frequency 
u>. I have assumed in the Fourier transforms throughout 
the previous sections a real frequency w. However, accord- 
ing to Landau's rule this must be understood as a lim- 
iting value of frequencies with negative imaginary parts 
A. The reasoning is that only perturbations, which grow 
in time from minus infinity are physically relevant, and, 
if perturbations with other imaginary parts are consid- 
ered, these must be analytical continuations of the former 
(cf. Binney & Tremaine 1987). According to the definition 
of the Fourier transform used here the functional depen- 
dence of the perturbations is e wt cx e~ At , so that negative 
imaginary parts A indicate growing perturbations. This 
implies that in the integral (51) the integration contour 
must be chosen in that way that it encompasses the poles 
of the integrand in the complex domain of r\ below the real 
axis. With these caveats in mind the integral (51) can be 
evaluated as 



drj e 



/ iA( n ' T n) V Tn 



V' 



/+ oc +7 

= 2ite^ iAn 8 (A) + ne^ iAn lim / drj'— - 



= 2ne^ lAn 5 (A) =p ne^ iAn 

„ [ + °° , , sin (At/) 
iir + iV / drj' 

J — 00 V 



(52) 



where rj' has been substituted for rj' = rj ± n. V indicates 
the Cauchy principal value of the integral, and if the latter 
is evaluated as 7r«sign(A) (Gradshteyn & Ryzhik 1965, 
equation 3.722), the integral (51) is equal to 



27re =F<An (J(A)T*n«(A)) , 

where u (A) denotes the unit step function, 



u(A) 



Inserting this back into equation (50) gives 



B n ( iui 




(53) 



(54) 



A k 

5(A) cos(nA) 




— n sin(nA)u(A) 
+zncos(nA)u(A) 



(55) 



where the upper row refers to the plus sign and the lower 
to the minus sign in equation (50), respectively. The case 
of positive wave numbers k' y can be treated in an analogous 
way, leading to a sign reversal, 



B K 



iuj 



(5(A)cos(nA) -nsin(nA)u(A) 
+zncos(nA)u(A) 



(56) 
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In order to assemble the results I consider first terms and 
proportional S (fc x — k x ) or (5(A) and integrate immedi- 
ately over k x , which contracts the integrand to fc x = k' x . 
Equations (48), (49), and (55) together with (47) give 



4irG k S 

OO 

■{I (A (k x = *£))- 2 £ (-l) n In (-A (fc x = k'J)} 



$ k < [2tt (-2B) - 2tt (-25) e - A (^= fe i) 



Since 



OO OO 

£ (-if In (-A) = £ Jn(A) - - (e A - Jo (A)) 



(57) 



(58) 



n=l 



(Abramowitz & Stegun 1968, section 9.6) , the term (57) 
cancels out exactly. This leaves according to equations 
(47) and (56) 

4ttG k So r+ °° 



.w, > / '^.'I'k 

2tt -2B k< J.^ 



•exp 1 



2 + + 4B 2 k ' 



B K 



■ exp — 



(59) 



•2 [ ]T(-n) sin (»A) cos (n («; - u;')) (-l) n J„ (-A) 

71=1 

-ncos(nA)sin(n(w-u7)) (-l) n / n (-A) u (A) . 
The series 

OO 

^nsin(n(A + w-uJ'))/ n (A) (60) 

n=l 

can be simplified by using the identity 
n/ n (A) = | (/»_! (A) - J n+1 (A)) (61) 
and rearranging it as 
^io(A) sin(A + w-uJ') 



Acos(A + tU-tU') = -| 



(64) 



{Mi 



7777 / '"" ) 



— ^(^-fc x )sinA} 



so that the multiple integral (38) can be put into the final 
form 



2ttG£o 
' An\k'J 



dk x <S> k {lk x k' x + j^k' y ) sin 



k [k x fc x ) 
2Ak' y 



+ ^ fc y( fc x- fc x)cOS 



2B 
• exp 2 



2Afc' 



}exp- 



. /c /c x 

-2<x> — 

2^ 

x 4B 2 y 



CQS K ; fee) + _*_ fc , (fc , _ fcx) gin « (fc x fcx) 



2Afc' 



2Afc', 



.(65) 



In expression (65) positive wave numbers k' y are assumed. 
In the case of negative wave numbers k' y the integration 

limits are reversed to dk x . 

The fundamental integral equation is thus a Volterra 
equation, 



#k' = / X dfc x /C(A; x ,fc x )$ k 

J — OO 

with the kernel 



(66) 



K, {k x , k x ) 



{ M x 



ttGSo 



AK l fc yl ^l + k^ 



exp — «cj- 



2Afc; 



—A:' 2 
4B 2 y 



K (^x ^"x) 

sin 

2AM, 



• exp 1 



2 r 1 

- ( A: 2 + k 



X J + 4B 2 



M x + -^K e2 



■ COS 



K (fc x fc x ) 

2AkJ y 



2B 



K (K - fc x) sin 



k (k' x k x ) 



2AM, 



■ (67) 



The integral equation (66) is equivalent to that obtained 
by Julian & Toomre (1966) although it is derived and 
presented here in a different way. 



+A sin (A + w - w) ^ I n (A) cos (n(A + w- w')) 



71=1 



A 



7. Discussion of the Volterra equation 

Surprisingly the Volterra equation (66) has no solution. 
= — sin (A + w — w') exp (A cos (A + w — w')) (62) This can be seen by a Fredholm discretization of the equa- 
tion (Morse & Feshbach 1953). The integral equation is 
(Abramowitz & Stegun 1968, section 9.6). Furthermore, then replaced by a set of algebraic equations with a trian- 
using the definitions of A, w, and w' I find the relations gular coefficient matrix. The kernel (67) vanishes on the 

diagonal, k' x = k x , so that the diagonal of the coefficient 



A sin ( A + ZJ - ID 1 ) = -| 



matrix is given by the identity matrix from the left hand 
side of equation (66). The determinant of a triangular ma- 
tt; , , ( , k 2 ,2\ trix is equal to the product of the diagonal elements of 
2B ky ^ x ~ k ^ ° OS A + I fcxfcx + 4B 2 k y ) Sm A ^ ^ 63 ^ the matrix, which is equal to one in this case. This means 
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that the set of algebraic equations and thus the integral 
equation have no non-trivial solutions. According to the 
theorem of Frcdholm the inhomogenous Volterra equation 



hi 



dk x K, (k x , k' x ) $ k + r k i 



(68) 



has a unique solution, if the corresponding homogenous 
equation has no solution. The solution can be expressed 
with the resolvent operator 1Z as 

rK 

*k' = rv + / dk x U (k x , k' x ) r k , (69) 

J — OO 

where the resolvent operator is given by a Neumann series, 

oo 

K(k x ,k' x ) = Y / K {n) (kx,K) (70) 
with 



71=1 



/C (1) {k x ,k' x ) = ]C{k x ,k' x ) ,and 

/C (n) (k x , k' x ) = / dk'^lC (*£', k' x ) KS n -V (k x , k'i) . (71) 

The inhomogcncity of equation (68) represents pertur- 
bations of the shearing sheet, which have not been in- 
cluded in the analysis so far. As will be discussed in the 
next section, these may be due to external forces or mas- 
sive perturbations within the disk, or any initial non- 
equilibrium states of the disk. Equation (69) describes 
then the response of the disk to these perturbations. 
If they are absent, the shearing sheet does not develop 
any inhomogenous spatial structures and perturbations 
of the velocity field. There is a technical point concern- 
ing the functional dependence of the resolvent operator 
on the frequency uj, which will become important later 
on. Equation (67) shows that the kernel JC (k x , k x ) is pro- 

k' — k 

portional to exp — iui 2 x ^ fc/ x ■ This proportionality is main- 
tained under the iterations of equation (71), so that the 
resolvent operator has the same functional dependence as 
/C, K (k x , k' x ) oc exp -iw k *~^ . 

When deriving the Volterra equation (68), I have 
not considered the case k' y — 0, which corresponds in 
disks with polar geometry to ring-like perturbations. It is 
straightforward to show using the formalism of section (4) 
that in this case the disk response to a Fourier component 
of the general potential perturbation is the correspond- 
ing Fourier component of the general disk response, so 
that they form conjugate pairs for each wave number vec- 
tor k' = (k' x ,0). When self-gravity is taken into account, 
the result is a dispersion relation lo (k x ), which is formally 
identical with that of Lin, Yuan, & Shu (1969), although 
it is valid here only for k' y = 0. Thus the case of k' y — is 
not generic for cases of even small non-vanishing circum- 
ferential wave members. The same holds true for Oort's 
constant A. The dynamical behaviour of a rigidly rotat- 
ing shearing sheet (A = 0) or star disks, in general, is not 
generic for the behaviour of differentially rotating disks. 



8. Initial value problem 

If the shearing sheet is initially - either at time minus in- 
finity or at some t = to - not in equilibrium, this must 
be taken into account explicitely in the Fourier transform 
of equation (20) with respect to the time coordinate. The 
Laplace transform of the time derivative of the perturba- 
tion of the distribution function gives in this case 

dt 

= -5f (t ) e -^o 



dte 



Sf (to) e 



to 



dt i 



to 



iiv5f u 



(72) 



This means a further inhomogenous term in equation (23). 
Equation (23) is linear, and therefore the solution of the 
extended equation can be constructed by superposing the 
solution of equation (23) in its original form as derived in 
section (4) and a solution of 



dSfu, 



fl ro) 



B dw2 
-5f(to)e- wt « =0. 



(73) 



In this way both inhomogenous terms of the extended 
equation (23) are taken into account. First, the inhomoge- 
nous term of equation (73) is Fourier transformed 



dfc y (5/ k (to) 



S1I11D1 




ifc. 



V2^7T 
2B 



■ COS W\ 



(74) 



Again the background distribution function / is split 
off all terms in equation (73), in particular 5f k (to) = 
/o/k (to), and a Fourier transform of equation (73) with 
respect to the angle variable W2 leads to 

r-\-oo 

rffc x /k (to) e- wta 



f 

j — ( 



• exp ik x 



J 2 - iloro 



-2B 



2Ji 



■ smioi 



• exp — ik y — — — cos wi = . 
2B 



(75) 



The derivation of the solution of equation (75) proceeds in 
exactly the same way as with equation (23) in section (4) , 
so that I report only the result using the notation from 
above 



f2.u 



i r +oc 

K J -oo 



■ exp i ( k 



f 



dk x f k (to) e~ wt " 



-2B 



dw[ exp i 



r/w 1 



2 sin (tttj) 

£ sin (w' x + wi — w) 



(76) 



B. Fuchs: Density waves in the shearing sheet 



The disk response to the potential perturbation is now where the w-dependent term has exactly the same form 

as in the kernel function JC. The general solution of the 
integral equation (68) has then the specific form 



given by 

Sfv = fo {fi,u> + f2,u>) cxp i (ojt + k y w 2 ) . 



(77) 



The second part of the disk response has also to be con- 
sidered in the Poisson equation. Again this can be done 
separately and the expression analogous to (38) is given 
by 



4ttG 



(2tt)' 
So 



dwi 



dw2 / dJ\ 



dJ 2 



-2B 



cxp ~ Ji 



dky exp(ifc y u>2) 



+oo 



dfex/k (*o) e- luJt " 



■ exp |^ ik x 
r+7T 



J 2 — ^0 r 

-2B 



— oo 

i 



2 sin (ttt]) 



i: 



■ cxp — I 



dw[ exp i [r)w[ — £ sin (w[ + Wi — w) 
i Ji - ^o^o , , / /2Ji 



-2B 



+ K 



sin wi 



■ expz 



, , ^/2kJ x 

fe y w 2 - fey 2B COSWi 



(78) 



The quadratures can be carried out in the same way as in 
section (6) and I report again only the result, 



2ttG£ 



• exp j 

K 



K 



rffc x /k x ,k' (to) e " exp -iu 



k x fc x 
2Ak' 



4B 2 



AB 2 



k x k x + ^ y ) C0S 



K (fc x fe x ) 

2AM, 



K (k'y- /Ox) 



2AM, 



(79) 



where positive wave numbers k' y are assumed. In the case 
of negative wave numbers k' y the integration limits have 
to be reversed to dk x . Thus the inhomogenous term of 
the integral equation (68) can be written as 



fV = 
with 



dk^C (fcx, k' x ) / kx ,k' (to) e 



-iujto 



(80) 



C (fc x , k x ) — 
• ex P-^2 



ttGSq 



k v 



2 



: exp — iw- 



2AM 



2 , «" ^ 



4B 2 



+ 



k v 



« 2 ,/2 

4_B 2 y 



( fc x fc * + ^2 fc 'y 



cos ■ 



K (Mx — ^x) 

2A/t: 



K (fc x k x ) 



2Ak' y 



(81) 



f- 

J — OO 



dkyrC (fc X) *x)/kx,k; (*o)e" 



(82) 



+ y X dk x TZ(k x ,k' x ) I ' dKC(k'Mf KK (t )e-^ ta . 



The potential (82) is still defined in Fourier space, but, 
when transformed back to time and spatial coordinates, 
gives the description of the evolution of the disturbance 
of the shearing sheet. Illustrative examples are described 
in the next section. 

It is instructive to consider the behaviour of the solu- 
tion (82) in the limit t — > oo (k' y > 0) , 



lim — 

t^oo 2tt 



lim 

t — >oo 



due^u'u 



*x/k x ,k' (*o) 



(83) 



■2Ak' y S (k x -k' x + 2Ak' y (t - to)) t (k x , k' x ) 



t—>OG 



/k' pk x 
dk x n (fc x , k'j / dk'lc, (k'l k x ) f KK (to) 
-oo J — oo 



■2Ak' y S (k£ ~k' x + 2Ak' y (t - to)) 
lim 2Ak' y £ (k' x - 2Ak' y (t - t ) , k' x ) 



■/ki-2AkJ r (t-t ),kJ, (h) 

r k '~ 

+ lim 2Afc' 

£— >oo y 
7ki-2Ak'(t-t ),k' (*o) 



& x 7e (fc x , fc x ) £ (fc x - 2Afc; (t - t ) , fcx) 



where operators with an overhead tilde denote again oper- 
ators with the w-dependent terms (cf. section 7 and equa- 
tion 81) split off. As can be seen from equation (81) the 
C functions decay as 

lim I (k' x - 2Ak' (t - to) , k' x ) 

= lim I (k' x - 2Ak' y (t - to) , fcx) 



a lim exp--^ (2Ak' y (t - t )) 2 = , 



(84) 



which means that the perturbation will die out eventu- 
ally. This was first noted by Julian & Toomre (1966). The 
physical reason is that the effective radial wave length of 
the perturbation, n/ (Ak' y (t — to)) , becomes much smaller 
than the typical epicycle size of the stellar orbits, mea- 
sured by ct u /k, so that the disk can no longer sustain 
such a perturbation. This is not the case for gas disks 
(Goldreich & Trcmainc 1978). 

9. Illustrative examples 

The characteristics of the solutions of the Volterra integral 
equation (68) become immediately clear, when the shear- 
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ing sheet is assumed to be perturbed initially by a single 
sinusoidal wave with a wave number k in (fcj" > 0) , 



sh (to = o) = f&s (k x - 4 n ) . 



(85) 



According to equation (80) this implies an inhomogeneity 
of the Voltcrra equation 



(86) 



and the solution of the Volterra equation is given accord- 
ing to equation (82) by 



+ 



/ dk x n(K,k'jfg a jC(k?,k x ) 



(87) 



Transforming this back to time and spatial coordinates 
leads to 



i r + °° 

' '* 2ir 



*k' ki 



2AM, 



(88) 



f k * ~ ~ ( k ln — k' 

+f§ J dk x n (fc x , *4) c (fc x n , fc x ) s [t ' 



2Ak\, 



and then 



S$(x,y,t) = — dk' x e^ x+k y "1 * k , , k ,» (t) 

2n J-oo y 



Ak 



y f^{£ (k?X + 2 Ak™i) 



tt y 

rk i ^+2Ak i °t 

+ / dk x K (fc x , fc x n + 2Afcj?t) £ (4 n , k x ) } 



■ exp i [(k x n + 2Ak™t) x + k™y] , 



(89) 



where the overhead tildes mean again that the in- 
dependent terms have been split off the operators C and 
1Z (cf . the remarks at the end of section 7) . An analogous 
expression can be derived for the surface density of the 
perturbation ST. Inserting (cf. equation 34), 



k 27rG k ' 

into equation (89) gives immediately 

1 Ak 



(90) 



ST (x,y,t) 



ft* ^{k™ + 2Ak™t) 2 + k» 



2-kG it -~v 

rk™+2Akft 

•{£(/c x n ,/c x n + 2Afcj n t) + / dk x 

J hi? 

■U (fc x n , k™ + 2Akft) I (fc x n , k x )} 
• exp i [(k x n + 2Ak™t) x + k™y] . 



(91) 



The spatial pattern described by equation (91) is still a 
sinusoidal wave, however with shearing wave crests. The 



wave crests, which are oriented perpendicular to the wave 
number vector, are defined by lines 



(fcj? + 2Ak i y a t) x 



k™y = const. . 



which tilt at the rate 
d ( Ay" 



dt \Ax 



= 2A. 



(92) 



(93) 



This corresponds exactly to the linear shear flow (11) of 
the shearing sheet. The amplitude function in equation 
(91) describes the growth and damping of the perturba- 
tion. As will be illustrated below, the amplitude function 
shows also some oscillatory behaviour, so that the solu- 
tions of the Volterra equation have also to some extent 
the characteristics of moving wave trains. For this complex 
phenomenon Toomre (1981) has coined the term 'swing- 
amplification'. 

Unfortunately the amplitude function in equation (91) 
cannot be calculated analytically. However it needs only 
very modest numerical work to solve the original Volterra 
equation by Frcdholm discretization numerically. The co- 
efficient matrix of the resulting set of algebraic equations 
is already in triangular form and extremely well behaved, 
because the kernel function (67) falls off in both the k x 
and k' x directions, respectively, tapered by gaussians. The 
Volterra equation can be written in dimcnsionless form, 
if the wave numbers are measured in units of the critical 
wave number, 



^crit — 



2irGT 



(94) 



and the radial velocity dispersion a u is expressed in terms 
of the Toomre stability parameter Q 



(95) 



3.36GS 



The numerical coefficient in the argument of the exponen- 
tial function of the kernel (67) is then given by 



„2 crit 



0.2856Q 2 



(96) 



Thus the kernel JC and the C function are parameterized by 
the values of A/^Iq and Q. The numerical scheme to solve 
the integral equation is briefly described in the appendix. 

In Fig. 1 the amplitude function of equation (91) is 
shown for the case of A/CIq = 0.5 and Q = 1.4. The ampli- 
fication factors are drawn as function of k x = k™ + 2Ak™t 
abscissae for increasing circumferential wave numbers k™. 
The sinusoidal waves have been incited with unit ampli- 
tude and with an initial radial wave number k™ = — 2 fc cr it. 
As can be seen from Fig. 1 maximum amplification is 
achieved, if the circumferential wave number is about 
0.5fc cr it- This is the exact equivalent of the 'X=2' cri- 
terion of Toomre (1981). The maximum of amplifica- 
tion is reached, when the effective radial wave number 
is k x s w 2fc cr it, which corresponds to fairly open trailing 
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'spiral' arms. The duration of the amplification event is 
given by 



Ucff 



t = 



2Akl: 



(97) 



and is thus of the order of an orbital period 2tt/V, . It 
can be also seen from Fig. 1 that, if the circumferential 
wave number k™ becomes very small, which corresponds 
to tightly wound 'spiral' arms, the duration of amplifica- 
tion is fairly long and approaches the asymptotic WKB 
limit (Lin, Yuan & Shu 1969). 




Fig. 1. Amplification of shearing density waves, de- 
notes the - constant - circumferential wave numbers in 
units of fc cr it- &x denotes the effective radial wave num- 
ber and is given by fc™ + 2Ak™t, so that it increases from 
left to right during a swing-amplification event. The ini- 
tial radial wave number is chosen as k™ — — 2fc cr ; t , which 
corresponds to initially leading arms. The parameters of 
the disk model are A/Q = 0.5 and Q = 1.4. 



While in Fig. 1 the arms were initially leading, in Fig. 2 
the amplification of initially trailing arms is shown. As 
can be seen from Fig. 2 there is hardly any amplification, 
which illustrates the well known result (Toomre 1981) that 
the swing-amplification mechanism requires leading input 
in order to amplify density waves. 




Fig. 2. Same as in Fig. 1, but the initial radial wave num- 
ber is chosen as = +lfe cr it- The parameters of the disk 
model are again A/fl = 0.5 and Q = 1.4 



10. Parameter study 

The solutions of the Voltcrra equation are parametrized by 
the values of A/Qq and Q. The peak of maximum amplifi- 
cation in wave number space shifts, when Oort's constant 
A is varied. I have calculated a grid of models with varying 
values of A/CIq keeping the Toomre parameter constant 
at Q = 1.4 and have determined for each model the cir- 
cumferential wave numbers for which peak amplification 
is reached. These can be fitted empirically by a relation of 
the form 



1.932 - 5.186 



A 



+ 4.704 



A_ 



(98) 



over the range 0.1 < (j^) < 0.5, which corresponds to 
rising rotation curves. This domain is particularly rele- 
vant, because in spiral galaxies the optically visible star 
disks often extend only over regions, where the rotation 
curves are still rising. Relation (98) is in perfect agree- 
ment with the result of Athanassoula (1984), who sudied 
the influence of the shape of the rotation curve on swing- 
amplification by a different approach. 




Fig. 3. Same as in Fig. 1, but the parameters of the disk 
model have been changed to A/Qq = 0.5 and Q = 1. 



The amplification factor itself depends strongly on the 
value of the Toomre stability parameter. In Figs. 3 and 4 
amplitude functions are shown which correspond to Q = 
1 and Q = 2, respectively. As can be seen from Fig. 3 
amplification becomes very high, when Q approaches the 
limit of Jeans instability Q — 1, while in the case Q = 2 
the disk is dynamically too hot to develop any structure. 

11. External perturbations 

When the shearing sheet is subjected to an external po- 
tential perturbation <f> cxt , the sheet will respond to it and 
the response generates its own gravitational potential <& mt . 
The overall gravitational potential is obviously given by 



$ = $ lnt + $ e 



(99) 



This has to be inserted into the Boltzmann equation, 
which leads to the determination of the distribution func- 
tion of the perturbed sheet and eventually to the right 
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(104) 



where the w-dependent term of the resolvent operator 1Z 
has been written separately (cf. section 7). Carrying out 
the integrations gives 



$ (x, y, t) = $ cxp{i [k x n x + kfy 1 ; 



Fig. 4. Same as in Fig. 1, but the parameters of the disk 
model have been changed to A/Qq = 0.5 and Q = 2. 



hand side of equation (35) . Into the left hand side of equa- 
tion (35), however, only the gravitational potential gener- 
ated by the perturbed sheet, <& mt = $ — $ cxt , has to be 
inserted. If the external potential perturbation is Fourier 
transformed with respect to time and spatial coordinates, 
its contribution to the left hand side of equation (35) can 
be calculated using equation (36) leading to an inhomo- 
geneity of the Volterra equation (68) 



• cxp 



kip 

1 

2T2 



t + 



2Ak\? 



(105) 



If the time constant is small, equation (105) can be ap- 
proximately written as 



$ (x, y, t) w $ exp{i [h£x + k™y] 

+$o exp i [(k x a + 2Ak y n t) x + k™y] 
■ti(k x n ,k x n + 2Ak y n t) . 



} 



(106) 



He' 



(100) 



11.1. Excitation of a density wave 



If the shearing sheet is perturbated by an impulse of the 
form 



$ oxt = $ exp{i [k x n x + k™y] - 



1 



}■ 



(101) 



The first term in equation (105) represents obviously the 
initial impulse, whereas the second term describes a swing 
amplified density wave exactly of the kind already found 
in the treatment of the initial value problem in section (9) 
(cf. equation 89). 

11.2. Response to a point mass 

The gravitational potential of a point mass resting at the 
origin of the shearing sheet is given by 



where the time constant T indicates the length of its du- 



Gm 



ration, the sheet responds with a swing amplified density $ = 

wave. The Fourier transform of equation (101) with re- 
spect to time and spatial coordinates leads to an inhomo- 
geneity of the Volterra equation (68) of the form 

TV = (2tt)^ <J> T<5 (k' x - 4") S (k' y - fcj") e~^ 2 T 2 . (102 ) 



Gm 



\/x 2 + y 2 



(107) 



Its Fourier transform with respect to spatial coordinates 
can be calculated using formulae (3.754) and (6.671) of 
Gradshteyn & Ryzhik (1965) giving 



The solution of the Volterra equation is given according ^oxt _ 2irGm 



to equation (69) by 



(108) 



$ k , = (2tt)* $ T5 (k' x - *£) 5 (k' y - 4 n ) e-*' 
+ / dk x K (k x , k'J (2tt) 5 $ T(5 (fc x - k£) 6 (k' y - k™) 

J —CO 

(103) 



so that the inhomogencity of the Volterra equation has 
the form 



1 , ,2rp2 



rk' = - 



47t 2 Gto 



Transforming this back to time and spatial coordinates 
leads to 



+ k'i 



(109) 



$ (x, y, t) = $ exp{i [k x n x + k^y] 



1 



} 



According to equation (69) the gravitational of the poten- 
tial of the perturbation is then 



+ 



_L f + °° 

2W-oo ' 



1 

'27T 



dk'e 



due 2Ak " 



k' 



*w = - 



4iT 2 Gm 



dk x TZ (k x ,k' x ) 

) 

(2n)H TS(K-k x n ) 



( K dk x n (fc x , k'j A ^ Gm j ( i 

J — oo 



1.2 j_ y 



(110) 
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and the corresponding surface density distribution (cf . sec- 
tion 9) 



k/ = 2itm8 (ui) 

J — c 



dk x n (fc x , k' K ) 



27TTO(5 (to) 

■\Jk%. + k' y 



(111) 



Transforming this back to time and spatial coordinates 
leads to 



to r°° r°° 

£ (x, y, t) = m5 (x) 6 (y) + — r / dk' x / dk! 

4tH J_ oa J_ oc 

f 

j — i 



dk* 



\Jk%. + k' y 



K (fc x , fc x ) A k '- x+k 'y y ^ . 



(112) 



The first term of equation (112) represents obviously the 
imposed point mass, whereas the second term describes a 
stationary polarization cloud induced by the point mass. 
Exactly the same polarization clouds have been found by 
Julian & Toomre (1966) or Toomre & Kalnajs (1991). 
Fig. 5 shows a contour diagram of the surface density 
distribution of such a polarization cloud calculated nu- 
merically by solving the Volterra equation with an inho- 
mogeneity according to equation (109) by Fredholm dis- 
cretization. I find the same surprisingly high surface den- 
sity of the polarization cloud as Julian & Toomre (1966). 
In the example shown in Fig. 5 the effective mass of the 
imposed perturber is increased by a factor of about ten 
due to the response of the disk. 

11.3. Passage of a point mass 

The dynamical effect of a satellite galaxy encountering a 
galactic disk on an unbound coplanar orbit can be sim- 
ulated in the shearing sheet by considering a point mass 
moving at constant speed vq along the y-axis. The exter- 
nal potential perturbation is then given by 



Gto 



\Jx 2 + (y-v Q t) 2 



(113) 



Using again the formulae (3.754) and (6.671) of 
Gradshteyn & Ryzhik (1965) the potential perturbation 
can be Fourier transformed as 



k' 



27rGm 



ik' vot 



(114) 



leading to an inhomogeneity of the Volterra equation of 
the form 



47r 2 Gm 



+ k'\ 



:d(u> + k' y v ) . 



(115) 



The solution of the Volterra equation has then an analo- 
gous form to equations (110) and (111). Transforming this 



.5 - 



- 



-.5 - 



-1 




-.5 .5 1 

x/'Vrit 

Fig. 5. Response of the shearing sheet to an imposed point 
mass illustrated by a contour diagram of the surface den- 
sity distribution. The inner peak represents the delta- 
function like point mass spread out according to the finite 
spatial resolution of 27r/fc max = A cr ; t /8 (FWHM). Contour 
levels at multiples of 25 m/A 2 rit are shown, where to de- 
notes the imposed point mass. Negative density perturba- 
tions are drawn as dashed lines. 



back to time and spatial coordinates leads to 

to f + °° f + °° 
£(x,y,i) = mS (x) 5 (y - v Q t) + — / dk' x / dk' y 

I 



K \/k' 
dfc x - 



VI _i_ U2 



^(fc x ,/fc x )e- 4 [S( fe x-4) + fe>0t] 



j\k' x x-k' y] 



(116) 



where the explicit form of the w-dependent term of the re- 
solvent operator 1Z has been used (cf. section 7) . Again the 
first term of equation (116) represents the imposed moving 
point mass, while the second term is now proportional to 



Jk' (y-v t) 



and represents thus a travelling polarization 



cloud which is stationary in the reference frame moving 
along with the point mass. 

11.4. Response to white noise 

Toomre (1990) has argued that much of the recurrent spi- 
ral activity seen by him in his numerical sliding grid sim- 
ulations of the dynamical evolution of a star disk or by 
Scllwood & Carlberg (1984) in their numerical simulations 
is due to swing amplified random fluctuations induced by 
particle noise. In order to illustrate this effect I have sub- 
jected also the shearing sheet to such potential perturba- 
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tions. I consider for this purpose the Fourier transform of 
the Volterra equation (68) with respect to time, 



*k',t (*) 



dkjt (fc x , k' x ) $ 



kx,k' t+- 



+ rv, t , (117) 



where use has been made of the explicit form of the in- 
dependent term of the kernel function and the convolution 
theorem of the Fourier transform of products of functions. 
This can be rewritten with the help of equation (90) in 
terms of the surface density 



-"k't 



k' 



dk K 



+E e J 



k'.t ' 



JC (fc x , k' x ) 

(118) 



where E k *\ denotes the surface density associated with 
the randomly fluctuating potential perturbations. 

In the case of white noise the autocorrelation function 
of the fluctuating surface density of the sheet is given by 



(S (x,i)E* (x',i)) 



J d 2 k J d 2 fc'(£ k (t),£ k , (t)) 



■expi [(k,x) - (k',x')] , 
where 

(E k (i) , ££, (t)) = |E k (t) \ 2 S ((k - k')/fccrit) 

so that 

(E (x,t)£* (x',i)) 

= *4 it y d 2 /c|E k (f) | 2 exp -i (k, x - x') 
= {2*fkl it \^{t)\H{x-x')8{y-y'), 



(119) 



(120) 



(121) 



because the amplitudes E k are uniformly distributed over 
wave number space. If the amplitude of the surface den- 
sity fluctuations is expressed as fraction of the background 
density of the unperturbed sheet, eE , one obtains 



(E(x,t) E* (x',t)) 



_2y2 



fc c 2 rit |E k (t)| 2 Afc x Afc y , (122) 



where the delta functions of equation (121) have been ap- 
proximated according to the finite resolution of the nu- 
merical solution of the Volterra equation (118). The inho- 
mogencity of the integral equation is then given by 



£ 2 £ 2 







k 2 



■ cos (<p) 



(123) 



"critA^xAfey _ 

defincd here as a real quantity, and where <p denotes a 
random angle distributed uniformly over the range [0, 2tt\. 
Finally it is advantageous to write the integral equation 
(118) and expression (123) in dimensionlcss form by di- 
viding them by the background density, which in wave 
number space is given by E (5 (k) s=a E / (Afc x A/c y ). Fig. 6 
shows a typical solution of equation (118). The parameters 
of the disk model are adopted again as A/Q = 0.5 and 



Q = 1.4. The time steps at which solutions of equation 
(118) can be calculated are given by multiples of 



At = 



Afc x 
2AM, 



(124) 



depending on the circumferential wave number k' y . Thus 
solutions arc obtained at different time intervals. However, 
within less than an orbital period a quasi-stationary equi- 
librium is reached, which is shown in the second panel 
of Fig. 6. Each spike represents the amplitude y^Ei^ 2 of 
the superposition of the randomly incited swing-amplified 
density waves, all travelling from left to right at a veloc- 
ity of fc Xje ff = 2Ak y . The contour plot in Fig. 6 shows 
clearly the peak of amplitudes around (fc x ,fc y ) w (1.7 , 
0.5) fccrit- This corresponds exactly to the amplification 
factors shown in Fig. 1. It was already noted in the nu- 
merical simulations by Sellwood & Carlberg (1984) that 
the constant spiral activity leads to dynamical heating of 
the disk, because the stars are randomly accelerated by 
the shearing density waves. This implies a rapid increase 
of the velocity dispersion of the stars and thus of the sta- 
bility parameter Q (cf. equation 95), which will eventually 
suppress any spiral activity, if the disk is not dynamically 
cooled. Disk heating by swing amplified density waves in 
the shearing sheet will be discussed in detail in a further 
paper of this series (Fuchs 2001, in preparation). 

12. Summary and conclusion 

The objective of this paper is to present a new description 
of the dynamics of a shearing sheet made of stars. It is 
shown again that disturbances of the disk evolve always 
into swing amplified density waves, i.e. spiral-arm like, 
shearing density waves, which grow and decay while their 
wave crests swing by from leading to trailing orientation. 
Specific examples are discussed how such swing amplifi- 
cation events are incited. Only in response to imposed 
massive perturbers the shearing sheet develops stationary 
perturbations. 
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Appendix 

Comparison of equations (89) and (87) shows that the am- 
plitude function in equation (89) can be calculated directly 
from the solution of the Volterra integral equation (68) with 
the inhomogeneity (86), if fc x takes the role of fc x n + 2Ak' y n t, 
(fc x > A;™), and K, and £ are replaced by K, and C, (u) = 0), 

rk' 



dfc x /C(fc x , fc x )$ kx , k ^ + £(fc x n , *£) . 



(125) 



Equation (125) is discretized as (fc x 

fki+l 

dfc x £(fc x ,fc x )$(fc x ,fc;) = 

'k, 

1 [ K.(h + 1, fcjO$(fe + 1, fey) + K.(ki, k-)$(ki, ky) Ak ( 126 ) 



leading to 



k',-1 

EAk 
~2~ 

fcj — — oo 



£ffc + 1, fej)$(fci + 1, fey) 



+£(fe, fcj)$(fe, fc y )J + £(fci n , fcj) . (127) 

According to the definition of the kernel /C(fcj, k'A = 0, so that 
the set of algebraic equations (127) can be resolved recursively: 

k- < k^ : $(fej,fe y ) =0, 

fej = : ^ ( fej i ) = £> ( 5 ^x ) j 

k- = *i n + 1 : $(fej, k' y ) = ic(k™ - 1, fci" + i)$(fci n - 1, fey) 

+ /C(fc x -\ fci" + l)$(fe- n , k'y) + K(kl, C + lWkl, k'y) 

- 1 Ah 

+£(fc x n +l,fe ln +l)$(fc x n + l,fcy) =■ 

+£(k?,k? + l) 

= £(fc x n , fci n + 1) + iC{k™, k in + l)$(fe x n , k' y )Ak, 
*(fej,fe y )=£(fe x n ,fej) 

fcf-1 

J 

+ V t(h, kj)$(ki, fe y )Afe . (128) 

In the numerical realization of the scheme I have adopted a 
stepwidth of Ak = 0.1 fc cr it. 



